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© Display system interpolators. 



© An image interpolator implements an interpolation function providing real time, continuous zoom capability r to 
„ image dfsp.ay system. Output image pixels are obtained by interposing the values «J^» r «^s ?£- 
of the 2x2 matrix of pixels surrounding the point on the input image. The disclosed arrangement employs a b.- 
t£ ntlpotL algorithm implemented in the form of cascaded one-dimensiona, interpolate crcu.ts 
Magnification control is established so that a unit increment of the zoom controller, such as a cursor on a tablet. 
Ss in i ^ constant increase in the degree of magnification. The coefficients required for the mtorpotabon are 
generated in real time avoiding the need for time consuming table look-ups. 
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The present invention relates to display system interpolators in or for display systems for processing 
image data for presentation on a display terminal, printer, or similar output device. 

(mage display systems have been developed to display digital images captured through optical 
scanning, video camera input, or via image sensors. The most basic image display systems simply 
reproduce the image stored in the digital image memory. More advanced systems allow for the transforma- 
tion of the stored image through colour enhancement, magnification (zoom), rotation, or other transforma- 
tion. 

The implementation of image magnification or zoom has been typically accomplished in the prior art 
through the us8 of pixel replication. In a system that uses pixel replication, an image to be magnified by a 
factor of 2 will be created by displaying two pixels of a given value or every one pixel of that value stored in 
the .mage storage. This form of magnification has the affect of creating what amounts to a single large pixel 
m place of each original pixel. At greater magnification factors this leads to considerable jaggedness in the 
>^ out,ine 01 the ima 9 e - Pixel replication is particularly undesirable in an image dispfay system which allows for 

- the dis P ,a y of mu,ti P ,e shad es of grey or colours. In these systems, the replication of individual pixels 

., creates a grainy picture of low quality. 

Theoretical techniques exist for interpolating the colour or intensity (grey shade) of the pixels of 
magnified images. The magnification process results in more than one pixel being displayed for every pixel 
of data stored. The addition of a pixel between a white pixel and a black pixel is better represented by a 
grey pixel than by the replication of either the black or white pixel. The implementation of these 
interpolation techniques, however, has not proved to be practical due to the considerable calculations 
required to interpolate the multitude of pixels in a display image. When implemented in software, 
interpolation results in unacceptably slow response times to magnification requests. 

Tabata et al.. "High Speed Image Scaling for Integrated Document Management" ACM Confe rence on 
Qffice '"formation Systems . Toronto, June 1984, discuss the use of high speed scalinglechniques relying 
on table lookup and shift operations. The suggested techniques speed the process but still require the 
building and referencing of a table. Tabata et al. interpolate based on a subdivision that has an interval 
related to the numerator of the magnification factor. This causes a loss of efficiency when fine resolution on 
the magnification factor is required. In addition, the Tabata et al. technique is based on the period of the 
interpolated sample deviation to input intersample distance, which tends to have a long period further 
30 decreasing efficiency. 

It is therefore an object of the present invention to provide an image display system having a 
magnification or zoom capability with real time interpolation of the pixels. 

Accordingly the present invention provides an interpolator in or for an image display system compris- 
ing: 

as image storage means for storing an image comprised of a plurality of picture elements - pixels - organised 
in two dimensions; 

a plurality of buffer means for storing one pixel from each position in a first dimension; 
means for transferring the stored image into the line buffers; 

coefficient generating means, responsive to an input magnification or zoom factor indication, for generating 
40 interpolation coefficients: 

filter means for interpolating output pixels from the line buffer in terms of the generated interpolation 
coefficients; and 

selection means for successively connecting predetermined number of the plurality of line buffers to the ' 
filter means. 

45 Such an arrangement does not require the use of table lookups and can be produced in a configuration 
that is thought to minimise the number of circuits and therefore the cost required for implementation. 
It can be used to display high quality magnified images of colour and grey scale stored images. 
The present invention will be described further by way of example with reference to embodiments 
thereof as illustrated in the accompanying drawings, in which:- 
50 R 9- 1 's a block diagram of an image display system incorporating an image processor; 

Fig. 2 is a diagram illustrating the relationship between the input image and interpolated pixels; 
Fig. 3 is a graph illustrating interpolation in one dimension according; 

Fig. 4 is a graph depicting the relationship between the values of TB { and the approximation TAj; 
Fig. 5 is a graph depicting the values of TBj and TAj after subtracting the constant slope; 
Fig. 6 is a table of derived factors; 
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Fig 7 is a block diagram of one form of image interpolator; 

Fig. 8 is a diagram depicting the relationship between input pixels and interpolated pixels in one and 

two dimensions; 

Fig. 9 is a block diagram of a two dimensional filter; 
5 Fig. to is a block diagram of a one dimensional filter; 

Fig. 1 1 is a block diagram of a coefficient generator; and 
Fig. 12 illustrates an alternative method of generating error term d m . 
One form of image interpolator according to the present invention operates in an image display system 
as shown in Fig. 1. The image display device is typically connected to a host computer system through a 
to host communications interface 10. A microprocessor 12 is provided to control terminal functions including 
host communications, data acquisition, managing an image storage 16, and providing parameters to an 
image processor 14. Image processor 14 is responsible for accessing image data stored in .mage storage 
16 and placing it in frame buffers 18 and 20 for display upon the monitor 22. The image display system 
typically uses a dual frame buffer system so that one frame buffer, for example, frame buffer 18. is used for 
,5 disolay generation while the second frame buffer, e.g. frame buffer 20, receives the image data from the 
image processor 14. The assignment of the frame buffers is reversed when the second frame bu.fer has 
[ M been filled Image processor 14 performs image manipulation functions such as translation, rotation, colour 

transformation, and magnification or zoom as described for the present invention. Image data .s stored in 
image store 16. typically a random access memory. Optionally, a graphics processor 24 may be included in 
20 the image display system to generate graphic objects based on graphics orders. Peripheral controller 26 .s 
provided to control input and output devices such as keyboard 28. a graphics tablet 30. or a v.deo input 

deV 'lma 3 ge interpolation is a function of image processor 14. Processing of the image without translation, 
magnification, or rotation involves essentially the mapping of each picture element, or p.xel in .mage store 
as 16 to a pixel in frame buffer 18 or 20. Magnification or zoom results in the addition of pixels in the frame 
buffer so that the resulting frame buffer image is larger than the image contained in image store 16 by a 
given magnification factor. For example, an image stored in image store 16 that has four p.xels by four 
pixels, when magnified by a factor of 2, would result in an image of 8 pixels by 8 pixels in the frame buffer. 
Thus the image processor must create an image with more pixels than exist in the original input. 
30 The problem becomes one of determining the intensity or colour of the newly inserted p.xels. This is 
accomplished by mapping the additional output pixels back to the input pixels, and interpolating the 
intensity or colour of the added pixels from the colour or intensity values of the existing input p.xels. For 

example. in Fig. 2 the original image was represented by pixels P(1.1). P(1.2) P(3.4). The magnified 

image creates additional pixels located as shown by U(1.1). U(1,2),... U(3.4). The value of the intensity o 
os pixel 0(1.1) is a function of the intensities of pixels P(1.1). P(1.2>. P(2.1). and P(2,2). The creation of a real 
time magnification or zoom requires that the interpolation of intensity be performed as rap.dly as poss.ble 
with hardware of a reasonable cost. 

Interpolation can be extended to colour images with the data from several components (e.g. red. green 
and blue, or hue. intensity, saturation) of the input image being processed in parallel by multiple 
40 interpolators and the results for each pixel combined and stored in a single word of the frame buffer. The 
interpolators could use the same set of parameters and share the same control logic and some of the input 

and output logic. i 

The specification of the magnification factor is accomplished by the .mage d.splay term.nal operator 
through the use of tablet 30 or other similar device. The display terminal operator will first .nd.cate that 
45 magnification or zoom is required. Then by movement of the tablet cursor or similar action the user will 
specify the magnification factor to be applied. In the disclosed arrangement, the system .s configured so 
that a movement of the tablet cursor a certain distance results in a given magnification regardless of the 
present magnification value. Additional movement will result in additional magnification by linearly mcreas- 

so Tdiscussion of the mathematical basis of this arrangement will now be provided. The disclosed 
arrangement of the present invention employs a known bi-linear interpolation algorithm. Although this 
description deals mainly with this particular algorithm, the concept may be extended to many other 
interpolation algorithms. The intention is to provide a continuous zoom capability, which means that the 
magnification factor can be any number in a permitted range and need not be limited to integers. 
55 The disclosed arrangement is based on two fundamental approximations in the interpolation process. 
First the location in an x,y plane of an interpolated point is approximated to the nearest point in a 
predefined, two dimensional grid overlaying the input image. This approximation implies an error in the 
intensity value of the interpolated point. This error can be made as small as required by defin.ng a 
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sufficiently dense grid. Second, the magnification factor, G, is restricted to be in the form of: 
G = £ ; A £ B and 
G S G mox 

and in which A and B are positive integers, A is fixed in the specific implementation, and G max defines the 
5 maximum magnification allowed by the implementation. A virtually continuous zoom can be obtained by 
defining A large enough to allow a sufficiently large set of possible G values. 

The disclosed arrangement uses a bi-linear interpolation algorithm. An output point, such as pixel U(1,1) 
of Fig. 2. is computed as a linear combination of the four nearest points. The bi-linear interpolation function 
is separable into x and y components. Thus, the algorithm can be implemented as a cascade of two one- 
io dimensional interpolations. Referring to Fig. 2. the value of the intensity at point M a" can be determined 
based on the intensities at P(1,1) and P(1,2). Then the intensity at point K b" can be determined as a 
function of the intensities at P(2 f 1) and P(2,2). Finally, the value for pixel U(1,1) can be determined as an 
interpolation in the y direction based upon the values at points a and b. 

The interpolation process will be explained with reference to Fig. 3 which represents one-dimensional 
ts interpolation with a magnification factor G = 2.4. Fig. 3 illustrates the interpolation between four input points 
at X n through X n0 . The corresponding intensities are l(x n ) through l( x . 3 ). A magnification factor of 2.4 
results in the creation of eight pixels in the frame buffer output image based on the four input pixels. 

The problem is then to find the intensity values of the points U m , U m *, t ... such that the ratio of the 
intersample distance between the output points to the intersample distance between the input points is 
20 constant and equal to the magnification factor. For bi-linear interpolation, it can be shown that any 
interpolated sample, for example. U m . lies on the straight line connecting l(x n ) and l(x nM ). Or in other words 
that: 

(1) U m = U(X m ) = (1-B m ) l(X n ) + B m i(X n+t ) 
where X n <X m <X ftM 
25 and, 



(2) 



30 



B = 
m 



X - X 
m n 



X - X 
m n 



X ^ - X 
n+1 n 



with T s being the input inter-sample distance. B m is the ratio of the distance between the interpolated point 
35 and the previous input point to the intersample distance, e.g. B m = 0.25 indicates that X m is 25% of the 
distance between X n and Xo*i< 

An approximation of the interpolated value is developed by first arbitrarily specifying a uniform 
subdivision of the input inter-sample distance X n + t - X n into T parts where T is an integer. Instead of 
computing the interpolated value U m at X m , which is difficult due to the requirement for higher arithmetic 
40 precision, the value of U(X' m ) will be calculated with X m being the point on the subdivision grid closest to 
X m . Then let V m = U(X' m ) which is the interpolated value at the approximated location x' m . 

As the original input sequence is bounded (assume it ranges between 0 and 255 as in usual images), 
the maximum error in taking V m for U m is also bounded. The maximum error is a function of how dense the 
grid is defined. For T = 256, the maximum error is equal or less than one half the quantisation step of the 
45 input sample. 

The previous equation (1) can be rewritten as 

(3) V m = (1-A m )l(X n ) + A m l(X n+1 ) 

where A m is an approximation of B m expressed as: 

(4) A m = t ROUND(T*B m ) 

so The function ROUNO(n) gives the closest integer to (n). 

This equation implies that the value of B m is known, however the disclosed arrangement uses a 
recursion formula to determine the value of A m as follows: 

(5) A m = I MOD(TVi + h + d m , T) 
where 

55 (5a) h = FLOOR(T/G) 

MOO(a.b) is a modulo b, and FLOOR(n) the integral part of (n). 

The correction term d m is made equal to 0 or 1 , such that: 
(5b) A m = f ROUND(TB m ) 
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This formulation means that the value of A m is obtained from the value of the previous coefficient A m ., 
by adding an integer h. This value is occasionally corrected by adding one (d m ) so that the selected gr.d 
point is the closest to the theoretical position. The problem is reduced to finding the appropriate sequence 
of values for d m . The example of Fig. 3 corresponds to a situation where T = 8 and G = 2.4 (which implies 
T'G = 3 333 and h = 3). Note that the abscissae of the samples V m , V m . 2 , V m . 3 . V m . 5 and v m . 6 are 
obtained by adding h to the corresponding abscissa of the previous sample. However the abscissae of 
V m . , and V m . 4 are obtained adding h + 1 . 

Fig. 4 graphically illustrates the relationship between the values T'A,- and T'Bj where i = m-1. m. m + 1. 
. T'Bi defines the location of a point with respect to the subdivision grid. In the example above. T"B m = 
3 333 indicating X m is between the third and fourth points of the grid. T*A m approximates TB m to the 
nearest subgrid point. I.e. in the example T'A™ = 3. The values of TB-, lie on a straight line with a slope of 
T/G The values of T*Aj are shown as the approximations of T'Bi lying on the subgrid. 

Fig 5 shows a line obtained after subtracting a line with a slope equal to h (in this example. 3) from the 
bne o' Fg 4. The errors between the values of T'B f and FA, can be expressed in units of the subdivision 
o,,.;, Fro-r> th.s figure it is apparent that the computation of the A m values can be done using an algorithm 
$.m.iar to Broserham's algorithm used for vector generation in graphics systems. A description of the 
B«'»nnam algorithm is presented, for example, in Fundamentals of Interactive Computer Graphics, by J. D. 
For, ar*3 a VarOam. published by Addison-Wesley. 1984. pp. 433-436. This type of algorithm computes 
t»x pc*is c« a o J'rote grid as close approximations to a given straight line. The slope of the line must be 
a fit- -:«m f-rrc-.x <<?ss than 1.0. 

A- '.uv;j ««. it is desirable to be able to have a linear increase in magnification in response to a 
spoof..- tm^r,: ci movement of the cursor on the tablet or similar input device. This makes it desirable to 
sc.--,., r*. m*yi»icaiion factor, G, as a function of an index L such that as a certain image is being 
d.sou.c.J - magnification factor G(L) and L is incremented, the resulting new image increases by a 
25 tactc* tna: is contact and independent of the value of L. In other words: 

(6t Gil •hGiU = Q . 

w.tn 0 coog a cywtant greater than one and representing the magnification increment associated with a 

unit ircfvas-; in tr»c ndex L. Equivalently: 
(7) G * Q"L 05LS Lmax 

fie asdosed arrangement, however, defines the magnification factor in a slightly different way to 
facilitate tne implementation in hardware. This expression can be converted into the form of equation (7) as 
will oe d scussea later. The magnification factor expression as used'in the disclosed arrangement is: 



20 



30 



N/2 



35 



(8) 



G = 



N/2 - K/w 



40 



45 



50 



in which N is the screen size in the dimension being considered, x or y, K as an index and w is a parameter 
fixed in the implementation. 

The meaning of the term K/w is related to the effect of magnification on the display. If K/w equals 0, 
then G = I indicating that there is no magnification and that the stored image will be presented exactly as 
it has been stored. Thus, for a screen with 1024 columns (N = 1024) the resulting display will also occupy 
1024 columns If K/w = 1. then G = 512/511. i.e. the image has been magnified so that 1022 of the input 
columns occupy 1024 screen columns on the output. If the centre of the image is assumed to remain in the 
same position, then one column at the left and one column at the right of the image have been discarded 
through this magnification. For K/w = 2. two columns at each side of the input image are discarded. Thus, 
the parameter 1/w specifies the fraction of the column that is discarded on each side of the input image 
when K is incremented by 1. i.e. K must be incremented by w to discard one column of the column at each 
side. The value selected for w determines the granularity of the magnification. A value of w in the order of 
magnitude of G m « provides smooth changes in magnification as discussed below. 
The magnification equation can be rewritten as 



55 D/2 

(9) G = 

D/2 - K 
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with 

D = w*N 

Since w and N are both fixed in this implementation, D will be a constant in that implementation. With 
this new formulation, the value of h, i.e. the integer slope (from equation 5a), can be expressed as: 

(10) h = FLOOR(T -- ££I ) = T - CE!L(2KT/D) 

where CEIL(x) is a function that returns the closest integer greater than or equal to the value of x. The slope 
of the line, Fig. 5, can therefore be expressed as: 

(11) S = CEIU2KT/D) - 2KTD 

(D/2T) CEIU2KT/D) - K 



D/2T 

If D/2T is an integer, the above expression is a rational number and, therefore, it is suitable for a 
Bresenham's algorithm. Furthermore, if T and D are selected such that 
(12) D/2T = 2~r 
20 r being an integer, then 

<2**r) CEIL(K/2**r) - K 

(13) S = 



25 2**r 

and the implementation is simplified. 

As an example consider the case N = 1024, T = 256, w = 2 which results in D = 2048, r = 2 and 



30 



35 



4 CEIL(K/4) - K 



(14) 



The table of Fig. 6 shows the values of G. S, h, and the periodic pattern of d m as a function of K for this 
particular example. Note that the pattern of d m has a period with maximum length equal to 2"r. 

The initial value for d. i.e. dO, is a function of how close the actual coefficient AO is to the theoretical 
coefficient BO. Its value will be given later. 
40 As stated above, this specification of magnification factor G, can be converted into a form that is a 
function of an index L such that a unit increase in L results in a constant increase in the magnification 
factor. The conversion between the two forms is accomplished as follows. 

The function G(K) increments most rapidly for G = G max . For this value of G, the corresponding value of 

Kis 



45 



50 



(15) K 

max = (1 - 



G 

max 



Let 

(16) G = G^-l) 
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D/2 



D/2 - (K - 1) 
max 



be the allowed value of G that is closest to G mn . Defining Q as: 



G 

max 

(17) Q = 

G ' 



then 




09) K = 2 (1 - Q"-L) which gives K as a function of L. (The obtained value of K should be rounded 
because K must be an integer). 
25 The range of L corresponding to the range 0-K„, aK of K is 0 to L max with: 



LOG(G ) 
max 



(20) L 

max 



30 LOG(Q) 

L can be entered by the terminal operator using a tablet for example. The computation of K as a function of 
L must be done once per frame or about 10 times per second. The equation can be computed by the 

as system control processor, directly or via a look-up table. It is still possible to use a look-up table for the 
conversion of certain values of K only and to compute the values in between by linear interpolation <p.ece- 
wise linear approximation). TU - i ~t n 

For N = 1024, w = 32. and G ma . =32. is 1777 (1776.189) and Q = 1.0019S3. Th,s value of Q 
guarantees that, in the range from G = 1 to G = 32. going from L to L + 1 results in one column being 

40 discarded at both sides of the previous image. This granularity is assured by the value w = 32 and causes 
the magnified image to change smoothly as the magnification is changed. Smaller values of w would resul 
in a more coarse implementation. A change in magnification would cause a noticeable jump in the output 
image when several columns are dropped. A one-dimensional magnification in the X ax.s has been 
discussed but the same techniques can be applied with respect to the Y axis. 

45 The implementation of such an interpolator is illustrated in Fig. 7. The interpolator block diagram in Fig. 
7 is a component of the image processor 14 shown in Fig. 1. The interpolator receives data from .mage 
store 16 on input line 15. Data from image store 16 is latched in an input register 50 and is then stored in 
line buffers 1 through 3 under the control of switch 52. The image processor accesses data from .mage 
store 16 on a row-by-row or line-by-line basis. The pixels from the image store are read into input register 

so 50 eg. 4 pixels at a time. Switch 52 causes the data to be transferred into line buffer 1 54 into , consecutive 
locations. Line 2 data is transferred to line buffer 2 56 and line 3 to line buffer 3 58. The 2D filter 64 uses 
the data of the first two lines of the input image to compute as many lines of the output .mage as required 
by the magnification factor G in effect.The result is transferred to one of the two frame buffers as .t is being 
computed The fourth line of input received from the input store 16 is transferred to line buffer 1 54 after the 

55 computation using the first two lines of data has been completed. During this transfer, the contents of line 
buffer 2 56 and line buffer 3 58 are used by the 2D filter 64 to generate additional lines of the output .mage. 
The process is repeated in the same fashion until the entire output image is computed and transferred to 
the frame buffer. 
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The line buffer select is controlled by line 68 output from the coefficient generator. Address generator 
60 raises the data from the input image stored in the line buffers to be transferred to the interpolator. The 
address is incremented in response to a signal on line 67. 

As discussed above, the two-dimensional filter required to interpolate the data can be implemented in 
s two one-dimensional steps. Fig. 8 illustrates how pixel V(n,m) is computed as a function of pixels P(i.j). P- 
+ 1 ). P(i + 1 j). and P(i + 1 j + 1 ). In a first step pixel Q(i,m) is computed as a function of P(i.j) and P(i.j + 1 ). 
Pixel 0<i + i,m) is also computed as function of P(i + ljJ and P(i + 1.j + 1). In a second step pixel V(n.m) is 
computed from the computed values of Q(i.m) and Q(i + l.m). 

The 2D filter 64 is shown in greater detail in Fig. 9. Two ID filters 70 and 72 compute the samples Q- 
;o (i.m) and Q(i + 1 .m). A third 1 D filter 74 computes the output sample V(n.m). 

In an alternate embodiment. 1D filter 70 can be eliminated since the sample for Q(i.m) has been 
previously computed by the ID filter 72 when processing the previous input line and. therefore could be 
stored m a buffer containing the complete output line. 

Fig 10 is a block diagram that illustrates the implementation of 1D filters 70. 72 and 74 in the disclosed 
75 arrangement r.-ns block diagram implements the equation for V m such that: 
(21 > v„ = (l-A_)ljXn) ♦ Aml(Xn + 1) 
which on r«? Afin-jn as 
(22) V, = kXni • Am(i(Xn + 1H(Xn)] 
which focmr i multiplication per output sample. 
20 Assume; tr^: in* .cput sequence is in the range 0 to 255, that T = 256, and Am is represented in 9 bits 
(8 bits lot cc<?fi«>;rt. one additional bit to allow rAm = 256). the range of possible products is in the 
closed mvrvaj <6S?30 €0^80) which requires 17 bits to be represented. The output adder takes the 9 most 
significant uts of cocuct and the input sample (with an additional positive sign bit) to compute the 
output same* rho lOtr most significant bit of the product enters the carry-in input to round the result. The 
25 result is contained m t?* 8 least significant bits of the adder output. 
The errcr a: tro output is a contribution of two terms: 

1. Ce-:nout*on term irom rounding the result. In this example this is equal or less than 0.5 of the 
quantisation level 

2. Contributor term from coefficient approximation. This term is the product of the error in the 
30 coefficient times the difference between the values of the two input samples. The worst case for a given 

approximation error n Am occurs when one of the samples is zero and the other is 255. This term is equal 
or less than 0.5 quantisation level. 

The total maximum error in each dimension is one quantisation level. The total error in both dimensions 
is equal or less than :wc quantisation levels for the example being considered. 
35 2D filters 64 and hence the 1D filters 70. 72 and 74 require the input of x and y coefficients to perform 
the interpolation. Thus there is provided a means of generating these coefficients in real time without the 
requirement for a look-up table. Coefficient generator 66 (Fig. 7) is illustrated in greater detail in Fig. 11. 

The value of Jhe coefficient Am can be obtained from the position of the computed sample in the 
subdivision grid (X m ). For the bi-linear algorithm, the value of the coefficient is: 



40 



- X (X' - X ) 

in n ran 

(23) A = 

m 



Vl " X n 



with T s the input inter-sample distance. The number T'A m gives the position of the sample with respect to 
X n in terms of units of the subdivision grid. Note that knowing T"A m allows computation of the coefficient 
so when other interpolating techniques are used provided the impulse response of the corresponding filter is 
known. For example, if a magnification by pixel replication is desired, the required coefficient can be 
computed as: 

(24) A m = SIGN(0.5 - A m ) 

SIGN (a) is zero if A is negative, and one if A is positive or zero. 
55 The computation of subsequent TA^ requires the computation of d m which can be obtained with a 
Bresenham's algorithm. Assuming that D/2T = 2~r, i.e. a power of two, the parameters of the BresenhanVs 
algorithm are: DX = 2""r 
DY = 2-r CEIL(K/2~r) - K 
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DY 5 ! D X ^s^^t^S of the line of Fig. 5. 00 is the initial error term of the Bresenham's 
algorithm, and Eo is the error in the position of the first sample in terms of units of the subd.v.s.on grid. i.e. 

X - X' 
m m 

(26) E Q = T ; (m = 0) 

X n + 1 " X n 

°' 5 The° expression for DO has been obtained from the equation that relates the error term, D, and the 
approximation error. E, in the Bresenham's algorithm: (27) D, = DX (2E r 1) ♦ 2DY 

One coefficient generator, as shown in Fig. 11. is required for each d.mens.on. T ,s assumed 1 to i be a 
power of two. for example T = 256. The system control processor 12 provides the parameters T AO h. DO 
and DY that apply for the whole frame. These parameters are stored in corresponds registers shown ,n 
£ 11 Component 108 represents a combinatoria! .ogic that provides a hard-wired left shift by one 
position, and EXORs the most significant bit with the output from INV ,14 (inverter, The CL 108 prov.des 
the parameters 2DY or 2DY-2DX to the adder according to the sign of the accumulator term <JATCH» 1UL 
The accumulator is initialised with DO at the beginning of each row or column^ Note that. as DX-2 r the 
term 2DY-2DX is equal to 2DY with the most significant bit changed. The EXOR function .s used for th.s 
purpose The initial address TA0 is loaded in the address accumulator ,04 at the beginning of each row o 
column The parameter h is added to the accumulator in every cycle by adder 100. The adder ^carry-in ,nput 
is fed with I from the accumulator 1,2. "Coefficient Logic" 106 is responses lor converting from the 
inerateTaddress to the actual coefficient value. It is not required for bi-.inear interpolation ,f T ,s a power 
InZ The adder ,00 cout (carry-out) signal is used to determine when a new input sample .s 
the carry-out of the x coefficient generator appearing on line 67 in Fig. 7 and the carry-out of the y 

C °tSSS:: "generate »e sequence d m which is practical foreman jJjjmJ . 
The r least-significant bits of K and a parameter related to the initial value of d m are stored m ^registe s K 
and D respectively. The counter is initialised with D at the beginning of a row or column *d .t ,s 
incremented as output samples are generated. The output d m is obtained as a function of K and the counter 
value using combinatorial logic 1 20. . . 

Assuming a 1024 x 1024 output image and that the interpolator ,s operating 75% of the time, a p.xel 
must be imputed every 75 nanoseconds, which is feasible using commercially availaWe components to 
cTstruct a cLit according to the disclosed arrangement. The remaining 25% of the time ,s reserved for 
possl bus contention of data coming from the image store, memory refresh, interpolator set up by the 
system control microprocessor or other functions. 
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Claims 

1 An interpolator in or for an image display system comprising: 
image storage means for storing an image comprised of a plurality of picture elements - p.xels - organ.sed 
in two dimensions; 

a plurality of buffer means for storing one pixel from each position .n a first dimension; 

f nr transferrina the stored image into the line buffers; 
ZZSt gSS^eans. response to an input magnificafion or zoom factor indication, for generatmg 

S ° ire^" "ItSating output pixels from the line buffer in terms of the generated interpofation 

2r e ctt n lan d s for successively connecting predetermined number of the plurality of line buffers to the 

"'^^interpolator as claimed in claim 1. including selection means for consecutively connecting two 
the plura.ity oHine buffers to the filter means at a time and frame buffer means for stonng the output of the 
filter means to accumulate the interpolated image. 

3. An interpolator as claimed in claim 2. wherein the coefficient generating means compnses: 
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logic means for determining the accumulated error associated with an interpolation coefficient; 

logic means for determining a new interpolation coefficient based on the initial interpolation coefficient and 

the integer interpolation increment incremented by unity, if the accumulated error is greater than or equal to 

zero; 

5 logic means for decrementing the accumulated error if the accumulated error is greater than or equal to 
zero. 

4. An interpolator as claimed in claim 2. wherein the coefficient generating means comprises: 
register means for storing values of 

an initial accumulated error. 
to an error correction term, 

an integer interpolation increment, and 
an initial interpolation coefficient; 

logic means for determining the accumulated error associated with an interpolation coefficient; 
logic means for determining a new interpolation coefficient based on the initial interpolation coefficient and 
is the integer interpolation increment incremented by unity, if the accumulated error is greater than or equal to 
zero; and 

logic means for decrementing the accumulated error if the accumulated error is greater than or equal to 
zero. 

5. An interpolator as claimed in any of claims 2 to 4, wherein the interpolation coefficients comprise: 

20 a first interpolation coefficient for interpolating in a first dimension; and a second interpolation coefficient for 
interpolating in a second dimension; and 
wherein the filter means comprises: 

a plurality of second filter means for interpolating between two pixels from one of the buffer means in 
response to one of the interpolation coefficients. 
25 6. An interpolator as claimed in claim 5. wherein the plurality of second filter means are connected such 
that 

a first of the second filter means is connected to a first line buffer for interpolating in a first dimension; 

a second of the second filter means is connected to a second line buffer for interpolating in the first 

dimension; and 

30 a third of the second filter means is connected to the output of the first and second of the second filter 
means for interpolating in a second direction. 

7. An interpolator as claimed in claim 6, wherein each of the plurality of second filter means comprises: 
difference means for determining the difference between the values of a first and a second pixel in the line 
buffer; 

35 multiplication means for determining the product of the difference and the interpolation coefficient; and 
addition means for adding the product and the value of the second pixel. 

8. An image display system comprising: 

an interpolator as claimed in any preceding claim; 

communication means for communicating with a remote data processing system; 
40 processor means for controlling the communication means and for storing an image into the image storage 
means; and 

data entry means for entering an interpolation command and a magnification factor; 

the coefficient generating means being responsive to both the interpolation command and the magnification 
factors. 
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® Display system Interpolators. 

0 An image interpolator implements an interpola- 
tion function providing rea" time, continuous zoom 
capability to an image d«sphy system. Output image 
pixels are obtained by inter ooiating the values of the 
colour or intensity of the 2*2 matrix of pixels sur- 
rounding the point on the input image. The disclosed 
arrangement employs a bi-linear interpolation algo- 
rithm implemented in me form of cascaded one- 
dimensional interpolation circuits. Magnification con- 
trol is established so that a unit increment of the 
zoom controller, such as a cursor on a tablet, results 
in a constant increase in the degree of magnification. 
The coefficients required for the interpolation are 
generated in real time avoiding the need for time 
consuming table look-ups. 
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